NAKED SINGULARITIES FORMATION IN THE GRAVITATIONAL 
COLLAPSE OF BAROTROPIC SPHERICAL FLUIDS 



ROBERTO GIAMBO, FABIO GIANNONI, GIULIO MAGLI, AND PAOLO PICCIONE 



O 
O 



Abstract. The gravitational collapse of spherical, barotropic perfect fluids is 
analyzed here. For the first time, the final state of these systems is studied with- 
out resorting to simplifying assumptions - such as self-similarity - using a new 
approach based on non-linear o.d.e. techniques, and formation of naked singular- 
ities is shown to occur for solutions such that the mass function is analytic in a 



^vq , neighborhood of the spacetime singularity. 



1. Introduction 



The final state of gravitational collapse is an open problem of classical gravity. It 
is, in fact, commonly believed that a collapsing star that it is unable to radiate away - 
via e.g. supernova explosion - a sufficient amount of mass to fall below the neutron 
star limit, will certainly and inevitably form a black hole, so that the singularity 
corresponding to diverging values of energy and stresses will be safely hidden - at 
least to faraway observers - by an event horizon. However, this is nothing more than 
a conjecture - what Roger Penrose first called a "Cosmic Censorship" conjecture 
ll34ll - and has never been proved. Actually, it is easy to see that one just cannot 
$h ! prove the conjecture as a statement on the mathematical evolution of any collapsing 

system via Einstein field equations, because in this case what is conjectured is 
baldly false: it is indeed an easy exercise producing counterexamples using e.g. 
negative energy densities or "ad hoc" field configurations. Thus, to go beyond the 
conjecture what is needed is a set of hypotheses, possibly based on sound physical 
requirements, which would allow the proof of a mathematically rigorous theorem. 
However, what turned out to be the truth in the last twenty years of research is that 
such a theorem (and, in fact, even the hypotheses of the theorem) is/are extremely 
difficult to be stated (see e.g. J23J). 

In the meanwhile, many examples of spherically symmetric solutions exhibiting 
naked singularities and satisfying the principles of physical reasonableness have 
been discovered. 

Spherically symmetric naked singularities can be divided into two groups: those 
occurring in scalar fields models [6, 8 1 and those occurring in astrophysical sources 
modeled with continuous media, which are of exclusive interest here (see ED for 
a recent review). The first (shell focusing) examples of naked singularities where 
discovered in dust models, numerically by Eardley and Smarr [ 10 1 and analytically 
by Christodoulou [5|. Today, the gravitational collapse of dust is known in full 
details E51 . 

The dust models can, of course, be strongly criticized from the physical point 
of view. In fact, they have the obvious drawback that stresses are expected to 
develop during the collapse, possibly influencing its dynamics. In particular, such 
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models are an unsound description of astrophysical sources in the late stage of the 
collapse even if the latter does not form a singularity: one can, for instance, regard 
a white dwarf or a neutron star as being an extremely compact planet, composed by 
a solid crust and a liquid (super)fluid core: such objects are sustained by enormous 
amounts of (generally anisotropic) stresses. It is, therefore, urgent to understand 
models of gravitational collapse with stresses. 

Recently, several new results have been obtained in this direction by considering 
systems sustained by anisotropic stresses (see e.g. [IT51 ITH1 l2T)l \2 1 1). Besides the 
details of the physics of the collapse of such systems, the general pattern arising 
from all such examples is that existence of naked singularities persists in presence 
of stresses: actually, we have recently shown that the mechanism responsible for 
the formation or whatsoever of a naked singularity is the same in all such cases 

na. 

In spite of the aforementioned physical relevance of anisotropic systems, it is 
beyond any doubt of exceeding interest the case of isotropic stresses, i.e. the grav- 
itational collapse of perfect fluids. In fact, for instance, the perfect fluid model is 
(in part for historical reasons) the preferred model used in most approximations of 
stellar matter of astrophysical interest. Unfortunately, although local existence ad 
uniqueness for the solution of the Einstein field equations has been proved 1281 1351 . 
very few sound analytical models of gravitational collapse of perfect fluids are 
known and, as a consequence, the problem of the final state of gravitational col- 
lapse of perfect fluids in General Relativity is still essentially open. Exceptions are 
the solutions describing shear-free fluids (see e.g. Il26ll27l0 and those obtained by 
matching of shock waves ll37ll : in both cases, however, the collapse is synchronous 
(i.e. the singularity is of the Friedmann-Robertson-Walker type) and therefore such 
solutions say little about Cosmic Censorship El l2"31 . 

There is a unique perfect fluid class of solutions which has been investigated in 
full details. This is the case of self-similar fluids, which has been treated by many 
authors since the pioneering work by Ori and Piran 031 (for a recent review see 
H). Self- similarity is compatible with the field equations if the equation of state 
is of the form p = ae (where p is the pressure, e the energy density, and a a con- 
stant). In this case the field equations reduce to ordinary differential equations and 
therefore can be analyzed with the powerful techniques of dynamical systems. Ori 
and Piran found that self similar perfect fluids generically form naked singularities; 
more precisely, they showed numerically that for any a in a certain range there are 
solutions with naked singularities. Recently, Harada added some numerical exam- 
ples which remove the similarity hypotheses ifTTl. 

These results clearly go in the direction of disproving any kind of censorship 
at least in spherical symmetry, since they show that naked singularities have to be 
expected in perfect fluids with physically sound equations of state. However, al- 
though being extremely relevant as a "laboratory", the self-similar ansatz is a over- 
simplifying assumption, and the general case of perfect fluid collapse remained 
untractable up today, essentially due to the lack of exact solutions. 

In the present paper we present the first (as far as we are aware) analytical study 
on the endstates of barotropic spherical fluids which circumvents this problem. To 
do this we use a combination of two new ingredients. The first is the fact that, in a 
suitable system of coordinates (the so-called area-radius coordinates) we are able to 
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reduce the field equations to a single, quasi linear, second order partial differential 
equation. As a consequence, the metric for a barotropic spherical fluid can be 
written, in full generality, in terms of only one unknown function. In this way 
the behavior of the null radial geodesies near the singular point can be analyzed 
in terms of the Taylor expansion of such a function. The second ingredient is a 
new framework for doing this analysis based on techniques for singular non linear 
ordinary differential equations lfT3lfi4ll . 

Our results here show the existence of naked singularities in barotropic perfect 
fluids solutions for which the mass function is analytic in a neighborhood of the 
center. 



2. Reduction of the field equations to a quasi-linear P.D.E. 

Consider a spherically symmetric perfect fluid. The general line element in co- 
moving coordinates can be written as 

(2.1) ds 2 = -e 2u dt 2 + e 2A dr 2 + R 2 (d6 2 + sin 2 #d0 2 ) 

where u, A and R are function of r and t (we shall use a dot and a prime to denote 
derivatives with respect to t and r respectively). Denoting by e and p the energy 
density and the isotropic pressure of the fluid, Einstein field equations can be writ- 
ten as 

(2.2a) = 47reR 2 R', 

(2.2b) * = -4tt P R 2 R, 

(2.2c) R' = Rv' + R'X, 

(2.2d) p' = -(e + p)i/, 

where \l/(r, t) is the Misner-Sharp mass, defined in such a way that the equation 
R = 2^ spans the boundary of the trapped region, i.e. the region in which outgoing 
null rays re-converge: 

(2.3) ¥(r,f) = - [1 - g^{d,R){d v R)} = | [l - (it")V 2A + {R) 2 e^ 
The curve t h {r) describing this boundary, i.e. the function defined implicitly by 

(2.4) R(r,t h (r)) = W(r,R(r,t h (r))), 

is called apparent horizon and will play a fundamental role in what follows. 

Initial data for the field equations can be assigned on any Cauchy surface (t = 0, 
say). Physically, the arbitrariness on the data refers to the initial distribution of 
energy density and the initial velocity profile, and is therefore described by two 
functions of r only. Data for R do not carry physical information and we parame- 
terize the initial surface in such a way that R(r, 0) = r. 

The data must be complemented with the information about the physical nature 
of the collapsing material. In the present paper we shall consider only barotropic 
perfect fluids, i.e. fluids for which the equation of state can be given in the stan- 
dard thermodynamical form: the pressure p equals minus the derivative w.r. to the 
specific volume v of the specific energy density e(v). We are going to work how- 
ever with the matter density p = 1/v and with the energy density e(p) = pe(l/p). 
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Therefore we are going to use in the sequel the equation of state of the fluid in the 
form (slightly less familiar than p = —de/dv): 

(2.5) p = pf- p - e 

Using the comoving description of the fields the matter density is proportional to 
the determinant of the 3 -metric, i.e. 



(2.6) p 



A-kER 2 



where E = E(r) is an arbitrary positive function. 

In order to simplify reading, we are going to develop in full details in the next 
sections the special - although physically very relevant - case of the linear equation 
of state 

(2.7) p = ae, 

where a is a constant parameter. However, in the final section, we will show how 
the results can be easily extended to (virtually) all the - physically valid - barotropic 
equations of state. 

In terms of the matter density eq. (12.51) implies e = p a+1 up to a multiplicative 
constant which however can be absorbed in the definition of E{r). For such fluids 
the field equation d2.2dl) integrates to 

(2.8) e u = p- a 

up to a multiplicative function of time only which can be taken equal to one by a 
reparameterization of t. 

We are now going to show that the remaining field equations simplify consider- 
ably (and actually the problem of the final state becomes tractable) if another sys- 
tem of coordinates, the area-radius ones, are used. The advantages of this system 
were first recognized by Ori ll32l . who used it to obtain the general exact solu- 
tion for charged dust. Subsequently, the area-radius framework has been success- 
fully applied to models of gravitational collapse and cosmic censorship (see e.g. 

Eimum). 

Area-radius coordinates are obtained using R in place of the comoving time. 
Denoting by subscripts derivatives w.r. to the new coordinates, we have \E' / = 
+ R' ^ jR , ^ = R ^ jR . Substituting in eqs (l2T2"ah . d2~2"ET) we obtain R' and p in 
terms of the mass: 

(2.9) R' = 

a + m R 



< 2 - io) '=b^r'' 

In writing the above formulae we have excluded the case a = 0. This case corre- 
sponds to the dust (Tolman-Bondi) solutions which is already very well known and 
will not be considered further in the present paper (see [25| and references therein). 
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Equation ( 12.31) can be used to express the velocity u = \Re~ v \ as 



u 



2^ 



(2.11) 

where we have introduced the function 
(2.12) Y = R'e~ x 

using d2~Sl> . (CT and dTTuT) we have 



Y 2 -l. 



(2.13) 



(a + l)p a 



This function plays the role of an "effective potential" for the motion of the shells. 
Notice that u is known when Y and ^ are; Y is known when E{r) is given and ^ 
is known. Thus, in particular, the initial velocity profile u(r, r) is known when the 
functions 



(2.14) 



\]> (r) = \[>(r,r) 



Yo(r) = Y(r, r) 



are known. It is therefore convenient to use Y as the second arbitrary function, 
eliminating E: 



(2.15) 



Y(r, R) 



ty, R {r,r)R 2 
*, fl (r, R) r 2 



Yn(r) 



where (12.101) and (12.131) have been used. 

We conclude that the metric for a barotropic perfect fluid in area-radius coordi- 
nates can be written in terms of the data and of the function \& and its first deriva- 
tives as follows: 



(2.16) ds 2 = 



dR 2 



R' 



2R'dRdr + [y \ (1 - ^)dr 2 



2tf. 



+ 



+ R 2 (d9 2 + sin 2 9d4> 2 ) 



where u, R' and Y are given by formulae (12. 1 Ob . (12.91) and (12.151) above. By 
a tedious but straightforward calculation the remaining field equation can be re- 
arranged as a second order equation for Remarkably enough, this equation is 
quasi-linear. In fact, the following holds true: 

Theorem 2.1. The Einstein field equations for a spherical barotropic fluid in the 
coordinate system (12.161) are equivalent to the following, second order PDE: 



(2.17) 



atf j . R R + 2&tf, P j l + c* rP = d, 
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where a, b, c, d are functions ofr, R, \l/ r , given by: 



(2.18a) 

(2.18b) 
(2.18c) 



1 



2 



'- [Y u)h 



1 — a 



a^ 2 



Y 



(a + l)V tR (Y 



u 



(2.18d) d=- 



a$ + R^ <R 




+ 



Y' 1 
Y a + 1 % (a 



2a 



IV 



Remark 2.2. Equation (12.171) must be supplemented with a set of data on the surface 
R — r. Since 

> 2 

(2.19) ac-b 2 



1 

a 



Y 



It 1 ® 



the character of the equation is determined by the sign of a. In particular, the equa- 
tion is hyperbolic for positive pressures and elliptic for the negative ones (recall 
that a = is excluded). For physical reasons, however, we consider here only 
the hyperbolic case (see next section). The initial data for equation (12.171) are thus 
given, in principle, by two functions. The value of ^ on the data corresponds to the 
physical freedom of assigning the initial mass distribution, while the first derivative 
can be calculated using eq. (12.91) evaluated on the data. On R = r one has R' — 1 
and therefore: 

(2.20) fl (r, r) = ^-^,,(r, r). 

a + 1 

Remark 2.3. A perfect fluid solution need not form a singularity: one can have 
oscillating, regular spheres as well. This poses the problem of characterizing the 
space of initial data w.r. to the final state (regular or singular). As far as we know 
this problem has never been studied (of course, it raises the issue of global existence 
that, as known, is extremely difficult) so that results like those known in the case of 
Einsten-Vlasov systems, for which 'small' (in a precise analytical sense) data lead 
to globally regular solutions [ 36 ] are not available here. In what follows, we are 
not going to address this problem. Therefore, we proceed further considering those 
data that lead to singularity formation with analytic mass function. It is, at present, 
unclear the degree of genericity of such data within the whole space of available 
data, and this will be the subject of future work. 



Remark 2.4. Equation (12.171) becomes degenerate at the sonic point, when the rela- 
tive velocity of the fluid equals the speed of sound. The behavior of the solutions at 
the sonic point is quite complicated, and not all the solutions can be extended. The 
problem of characterizing the structure of the space of the solutions is extremely 
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interesting. As far as the present authors are aware, such an analysis has been car- 
ried out in full details only in the self-similar case E1IS1IH1- I n me present paper, 
however, we are interested only in singularities which arise from the gravitational 
interaction. 



3. FORMATION AND NATURE OF SINGULARITIES 

3.1. Physical requirements. We are going to impose here strict requirements of 
physical reasonableness. First of all, we impose the dominant energy condition, 
namely, energy density must be positive and the modulus of the pressure cannot 
exceed the energy density (so that — 1 < a < 1). We consider, however, only the 
case of positive pressure. It must, in fact, be taken into account that, while tensions 
are common in anisotropic materials, a perfect fluid can hardly be considered as 
physical in presence of a negative isotropic pressure. 
Therefore, a > and ( 12.101) imply that 

(3-D V >R {r,R)<0, Vr>0, W?e[0,r], 

and since we want R' > to avoid shell-crossing singularities (see below), it must 
also be, from (12.91 ). 

(3.2) ty !r (r,R)>0, Vr>0, Vffe[0,r]. 

As mentioned above, we require the existence of a regular Cauchy surface (t = 0, 
say) carrying the initial data for the fields. This requirement is fundamental, since 
it assures that the singularities eventually forming will be a genuine outcome of 
the dynamics. It is easy to show that, with the equation of state used here, it is 
equivalent to require the matter density to be finite and non vanishing on the data. 
Due to eqs. (12.101 ) and ( 12.201 ) we get 

(3.3) lim *' r ^ ,r) £(Q,+oo). 

Since area-radius coordinates map the whole set {(t, 0) : t < t } into the point 
R — r — 0, one may ask whether this may give rise to some kind of contradiction, 
that is whether the hypersurface {R = r} fails to be regular. However, note that the 
coordinate change, restricted on the initial data hypersurface, is regular up to the 
centre, since the generic point (0, r) in comoving coordinates is mapped onto the 
point (r, r) in area-radius coordinates. Moreover, we are going to put analiticity of 
the data into play. In a neighboorood of the center, this property has to be checked 
using a cartesian system of coordinates, since even powers of r can give rise to loss 
of differentiability at finite order in such coordinates. 

To inspect this point we consider the whole set of Cauchy data for the fields. Let 
us choose a coordinate system on E in such a way that the embedding reads 



(3.4) 



E(a, 9, 4>) ^ M(r = a,R = a,9,(f>). 
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The induced metric and the extrinsic curvature (i.e. the second fundamental form) 
are respectively given by 

(3.5) dsl = . \ rda 2 + a 2 dn 2 , 

47TE{a)a z p{a, a) 

(3.6) = -^^r ( -±r) (a,a)da 2 -au(a,a)dQ 2 . 

It is now relatively easy to check that, if \^(r, R) is analytic and odd, and Y (r) = 
1 +0(r 2 ) is even, using (12.101) . (12.111 ) and (12.131) the above tensors on £ are analytic 
and even in r. 

This means that all the physical quantities give rise to analytic functions in carte- 
sian coordinates near the center. 

Finally, we require regularity of the metric at the center that is, in comoving 
coordinates: 

(3.7) R(0,t)=0, e A(0,t) = R'(0,t), 

for each t > up to the time of singularity formation t . 

The singularity forms whenever the denominator in (12.101) vanishes, that is R = 
0. This kind of singularity is called a shell-focusing singularity (we have excluded 
here, via equations (12.91 ) and (13.21 ) , the so called shell-crossing singularities at 
which the particle flow-lines intersect each other). In comoving coordinates (r, t), 
the locus of the zeroes of R(r,t) defines implicitly a singularity curve t s (r) via 
R(r, t s (r)) = 0. The quantity t s (r) represents the comoving time at which the shell 
labeled r becomes singular. The singularity forms if t s (r) is finite for each shell. 
In physically viable cases the curve t s (r) is strictly increasing and the center is the 
first point which can become singular. Regularity of the data then implies 

(3.8) lim t s (r) = t >0. 

In order to describe the singuarity formation at the shells r > by condition 
R = 0, from (12.101) we make the assumption 

(3.9) lrm — = -oo, 

R^0+ R 2 

for r sufficiently close to 0. Using the above requirements, toghether with (I2.12I) . 
we can also translate relations (13.71) in area-radius coordinates asking 

(3.10) lim Y(r,xr) = 1, Vi 6 (0, 1]. 

r— >0+ 



3.2. Taylor expansion of the mass. As said in Section 1231 in the present paper 
we assume analyticity of the mass function at (0, 0). It should be noticed that the 
'point' (0, 0) in mass-area coordinates 'contains ' both a regular part (it contains the 
data R = r as r goes to zero) and a part at which the spacetime becomes singular 
(as R goes to zero along the singularity curve, see next section). The mass function 
itself however satisfies an equation which is regular at the spacetime singularity, so 
that the assumption made here is exactly equivalent to that usually made on the data 
in other models of gravitational collapse. Such data can be taken to be analytic in 
cartesian coordinates near the center, as in 0, or simply Taylor-expandable up to 
the required order as in |25 1). In the present paper however we assume analiticity. 
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Moreover, coherently with our choice of initial data, we will assume odd-parity of 
the mass function. 

The following holds true: 

Proposition 3.1. The Taylor expansion of the mass function ty(r,R) has the fol- 
lowing structure 

(3.11) ^( r ,R) = ^(r 3 --^jR 3 ) + ^ + . . . . 

^ a ' i+j=3+k 

where k is an even integer, k > 2 and h is a positive constant. 

Proof. Odd parity of ^ and regularity condition (13.31) and (12.201) imply the Taylor 
expansion to start from third order terms. Therefore, one certainly has 

(3.12) *(r, R) = ^i/ Rj + ■■■■ 

i+j=3 

For the sake of convenience we now set, for each n > 0, 

(3.13) Mr)= lxV ^ J - B "^= E ^ijr j -\ 

i+j=3+n i+j=3+n 

so that the r n+2 's coefficients of Taylor expansions of ^, r {r, rr) and ^ ,#(?", rr) are 
A n {r) and B n (r) respectively. We recall that (13.31 ) implies A (l) > 0, and, from 
B (l) < follows. Using d2~T31) we get 



Y(r, rr) 



4>(i) 



5o(l)r 2 



at least for each r e (0, 1] such that B (r) ^ (but this polynomial can possibly 
vanish only for two values of r), and then (13.101) holds if 



(3.14) 
But 



Mr) 



B (iy 

Bo(r) 



Vr G (0,1] with B (r) ^ 0. 



B (l)r 2 2 *12 + 2^ 2 i + 3^3 
— r 



B (t) ^2i + 2^ 12 r + 3^o3r 2 ' 

and therefore if ^21 was not vanishing, the above quantity would tend to zero as 
r — > 0, which is in contradiction with (13.141) . Then ^21 = 0. A similar argument 
applies to ^12 to show that this quantity is zero as well. Finally, relation (12.201) 
imposes a constraint on A n (l) and B n (l): 

(3.15) -«A,(l) = (a + l)5 n (l), Vn>0. 

Using this equation for n = and setting h := 2A (1) we finally get formula 

<EEB. " □ 

Remark 3.2. A tedious but straightforward calculation shows that the Taylor ex- 
pansion (13.1 II) is compatible with (12.171) "in the Cauchy-Kowaleski sense" at any 
order, that is, the equation allows the iterative calculation of all the higher order 
terms once the data are chosen. Of course, we stress that this is not a proof of 
global existence up to singularity formation but only a - fundamental - consistency 
check for solutions here assumed a priori as regular. 
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Remark 3.3. The Taylor expansion given above excludes the self-similar solutions 
from what follows. It can, in fact, be easily shown that analyticity in self-similar 
variables leads to mass functions of the form ^ = rijj(R/r) where ^ is finite at 
R = 0. One recovers here a fact which is very well known in the case of dust 
spacetimes, where the self-similar mass profile has a constant \I/ (linear profile) 
while analyticity of the data for non self-similar solutions requires \& to start from 
cubic terms. 

3.3. The apparent horizon. A key role in the study of the nature of a singularity 
is played by the apparent horizon th(r) defined in (12.41) (see for instance |23 1). The 
apparent horizon is the boundary of the trapped surfaces, and therefore represents 
the comoving time at which the shell labeled r becomes trapped. In area-radius 
coordinates this boundary is defined by Rhor = 2\I/(r, Rhor)- Since ^^(O, 0) = 0, 
implicit function theorem ensures that the curve Rhor is defined in a right neigh- 
borhood of r = 0. In what follows, we shall need the behavior of this curve near 
r = 0. It is easy to check that R hor is strictly increasing and such that R hor {r) < r. 
Moreover it is Rh or (r) = 2\l>(r, 0), since from (12.41) it is 

R hor = 2* (r, 0) + 2R hor V >R (r, 0) + R 2 hor g(r, R hor ), 

where g is bounded and ^^(r, 0) is infinitesimal. Therefore, due to eq. (13.1 11) . we 
conclude that 

(3.16) R h (r) = hr 3 + . . . . 

Next section is devoted to the study of the nature of the central (R = r = 0) 
singularity. We restrict ourselves to this singularity since, in barotropic perfect 
fluid models with positive pressures, it is the only one that can be naked. This is 
easily seen using comoving coordinates. Indeed, a singularity cannot be naked if it 
occurs after the formation of the apparent horizon (i.e. it must be th(r) > t s (r)). A 
necessary condition for this is that the singularity must be massless (\I/(r, t s (r)) = 
0). But, due to equation (I2.2bl) . in presence of a positive pressure the mass is strictly 
increasing in a collapsing (R < 0) situation, while it is zero at the regular centre. 
The situation can be completely different if negative pressures are allowed: in this 
case non central singularities can be naked as well |9|. 

3.4. Nakedness of the central singularity. At the center (R — r — 0) the appar- 
ent horizon and the singularity form simultaneously and the necessary condition 
for nakedness is satisfied. The singularity will be (locally) naked if there exists a 
radial lightlike future pointing local solution R g (r) of the geodesic equation with 
initial condition R(0) = "travelling before the apparent horizon", that is - in area 
radius coordinates - R g (r) > Rhor( r ) for r > 0. We will study in full details only 
the existence of radial null geodesies emanating from the singularity. It can in fact 
be proved that, if a singularity is radially censored (that is, no radial null geodesies 
escape), then it is censored OTl H31 . 

The equation of radial null geodesies in the coordinate system (r, R) is easily 
found from (12.161) setting ds 2 = together with d# = d(p = 0: 

d-R a \P r / u 

(3 ' 17) IF^—itA'-Y 

Our main result can be stated as follows: 
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Theorem 3.4. For any choice of initial data Yq(t), ^o(r) for the Einstein field 
equations such that 

(1) the central singularity forms in a finite amount of comoving time, and 

(2) the Taylor expansion of the mass function is given by (13.1 II) . 

there exists solutions of (13.171) that extend back to the central singularity, which is 
therefore locally naked. 

To show the result we first need the following lemma. 

Lemma 3.5. Called t x {r) the curve defined by R(r,t x (r)) = xr 3 , there exists a 
x > h such that 



(3.18) 



lim t x (r) = t - 



Proof. It must be shown that for some x > h 

p a (r,a) 



(3.19) 



lim 



da = 0. 



With the variable change o = rr 3 the integral above becomes 

p a (r, r 3 r)r 3 ^ 2 ^ 



l (2#(r, r 3 r) + rr 3 (Y 2 (r, r 3 r) - 1)) 1/2 
and to prove (13.191 ) using Fatou's lemma it suffices to show that 



dr, 



(3.20) 



lim sup 

r->0+ 



p a (r, r 3 r)\/r 



\ 



2*(r,r 3 r) 



r^j + r Y 2 (r, r 3 r) 



1/2 



dr < +oo. 



/ 



We first notice that the quantity in square brackets at the denominator in the above 
expression must be positive for r small. This is to ensure dynamics near the central 
singularity (see, e.g., d2.HI) '). But, using (13.111) . it is 



2\P(r, r 3 r) 



(h-r) + 0(r 2 ), 



where 0(r 2 ) is infinitesimal uniformly in r (again, this notation means infinitesimal 
behaviour, uniform in r). Since r can be greater than h, then Y(r, r 3 r) cannot be 
infinitesimal as r goes to 0. Recalling Y = fj+^pc* > an( ^ exploiting (13.101) for x = 1, 
it is also a simple task to check that E(r) behaves like r~ 2 , 

E(r) V, r {r, r 3 r) = c + 0(r), 

and so p a (r, r 3 r) cannot be infinite as r approaches 0. The expression for p is given 
by (12.101) : for simplicity we compute p a+1 , using (13.1 II) : 



p a+ \r,r 3 r) = - Cl 



3 , a 
h- 



r 3 r) 



2 a + 1 t 2 



■ii 



cr 



+ ^ei + 0(r 



(^32 + 0(r)) + 0(r) 
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where c\ is the positive constant (A-kol) 1 . As said above, this cannot be infinite 
and therefore vanishes, giving for some constant c 2 



P a (r,r 3 r) = ^(b(r)+0(r))^, 

f a + l 



where b(r) is a regular function. This yields, passing to the limit r — > + , the 
following expression for the integral in (13.201) : 

(3-21) / 2 1 , V —^dr. 



2a 
f a + l 



(h — r)fe(r) "+ 1 + tCqc\t<*+ 1 



The term in square bracket at the denominator is bounded away from zero for r < h 
and so is for x greater than but sufficiently near to h. Recalling the bound a < 1, 
the above integral is therefore finite, and the lemma is proved. 

□ 

Remark 3.6. Let us observe that we have incidentally shown here that 

(3.22) *,ji(r, xr 3 ) = -a(x)r & + . . . , (l - y^j (r,xr 3 ) = d(x) + . . . , 

where a(x) and d(x) are some positive functions. 

Also observe that the same argument of the above lemma can be used to show 
that also t hor (r) tends to t as r — > + , that is the centre gets trapped at the same 
comoving time it becomes singular. 

Proof of theorem \n\ To show the existence of singular geodesies we use a simple 
technique developed earlier iTOl. First of all, we recall that a function yo(r) is called 
a subsolution (respectively supersolution) of an ordinary differential equation of the 
kind y' = f(r, y) if it satisfies y' < f(r, y ) (respectively >). Now, it can be shown 
lfl2ll that the apparent horizon Rh{t) is a supersolution of the geodesic equation 
(I3.17I) . The singularity is certainly naked if it is possible to find a subsolution R+(r) 
of the same equation which stays over the horizon. In fact, choose a point (ro, Ro) 
in the region S = {(r,R) : r > 0,R hor (r) < R < R+(r)}. At this point the 
(regular) Cauchy problem with datum R(r ) = R admits a unique local solution 
R g (r). Now the extension of this solution in the past cannot escape from S since 
either it would cross the supersolution from above or it would cross the subsolution 
from below. Thus it must extend back to the singularity with lim r ^ + Rg( r ) — 0. 

We now proceed to show that a subsolution always exist. For this aim, it suffices 
to consider a curve R x (r) = xr 3 , with x > h. Indeed, computing the righthand 
side of (13.171) for R x (r), using (13.221) . we get that R x (r) is certainly a subsolution 

of (EH) if 

(3.23) x < — ^— /\ 4 d(x), 

that is always satisfied, independently of x, for r sufficiently small. 

Therefore, if we consider the curve R x (r) for x > h sufficiently near to h, then 
Lemma l3~5l ensures that - re-translated in comoving coordinates - it emanates from 
the central singularity, and so the theorem is proved. 

□ 

We stress that the theorem holds for any solution satisfiyng (13.81) and (13.1 II) . 
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4. EXTENSION TO THE GENERAL BAROTROPIC CASE 

We are going to show in the present section that our main result, namely the 
existence of naked singularities, actually hold for the general (i.e. not necessarily 
linear) barotropic equation of state e = e(p) provided that a set of (physically 
motivated) requirements are satisfied by the state function: 

Assumption 4.1. We assume e = e(p) to be a C 1 function in [p, +oo) (where 
p > 0), such that e(p) > (= iff p = p). Recalling (H3l . that is p(p) = pj^ - e, 

we also assume p(p) is a strictly positive C 1 function with ^ > 0, except at most 
for a bounded interval [p, pi\, possibly coinciding with a single point, where p(p) 
can vanish. 

Remark 4.2. Observe that: 

(1) The assumptions made imply that 

(4.1) ^i(p) > ifp>p. 

dp 

and therefore e(p) is a strictly increasing positive function. 

(2) Differentiating (12.51 ) we have, where it makes sense, 

(A ™ d P & e 

(42) d-p =P dp-^ 

then e(p) is strictly convex for p sufficiently large, and so 

(4.3) lim ^(p) = +oo. 

p^+oo dp 

(3) The assumptions made imply the existence of linip^oo p(p). In addition, if 
the limit would be finite, say I, then we should have ^(p) < -(I + e(p)), 
and then e(p) < p + / by a simple comparison argument in o.d.e., which is 
in contradiction with (14.3b. Thus 



(4.4) lim p(p) = +oo. 

Remark 4.3. We stress that the above mentioned hypotheses are quite natural from 
the physical point of view. Besides obviously including the p = ae equation of state 
considered so far, they include, for instance, the equation of state of the perfect gas 
p(p) = K 2 p for which e(p) = K\p + K 2 p log p where K 1 and K 2 are positive 
constants (in this case one obviously has p\ = p = e~ Kl ^ K2 ). 

Einstein's equation (I2.2bl ) reads 



(4 - 5) p - toV 

Using it, together with (I2.2ab . (I2.5I) . (12.61) and the coordinate change formulae ^' = 
\l/ ir + R' ^f t n and ^ = R ^/^r, we obtain the general counterparts for equations (12.91) 
and (12.131) . namely 

*r P * r 



(4.6) R' 



m yR e + p AnR 2 {e + p)' 
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and 

£(r)tf r (r, R) 



(4.7) Y(r. 11) , 

(de/dp)(p(r, #)) 

Using these formulae, we can again express the metric in the form (I2.16I ). The 
crucial point is now that we can express ^(p) as a function of p(p) and, as a conse- 
quence, the dynamics of the system is, also in the general case, expressed in terms 
of the mass function and its derivatives only. To this end, consider the parameter- 
ized curve in 3? 2 



This curve is globally the graphic of a function ( = £(£), recalling that £ is a 
non-decreasing function of p (by the assumptions on p(p)), and that, by (14.2b . 
[dg/dp) — P 7^ 0. Since also ( is increasing by ^ of Remark l4~2l there exists 
lim p ^p+ C(p) = £ finite. Denoting by £ the finite number lim p ^p+ £(p), the func- 
tion may be prolonged up to the point (£ , Co)- Let also observe that this function is 
C 1 , for each £ > £ > where indeed is strictly positive by the assumptions made 
onp(p). 

Using this result, and recalling (12.5ft . one finds that Y in (14.71 ) (and then it in 
(12.1 lb ) can be expressed as functions of the data and of the mass function \I/(r, R) 
and its derivatives. Then, again, with some calculations one obtains a second order 
PDE that must be satisfied by \E'. As in the case treated so far, we consider only 
analytic solutions of this equation, and proceed to analyze the structure of the lower 
order terms of the mass profile. 

First of all, since R' = 1 on the data surface R = r, the expression for the initial 
energy is: 

¥, r (r,r) + ttji(r,r) 
<r,r) = — 2 • 

Imposing the regularity condition lim r ^ + e ( r ) r ) £ (0, +oo) and making reference 
to the notation used in Section I3T21 we get 

A (l)+5o(l) >0. 

Actually, Aq(1) ^ 0. Otherwise, e(r, r) + p(r,r) — > as r — > + , since e(r, r) + 
p(r, r) = ^f^P- = + o(l). But (e + p)(p) is a strictly increasing and non 
negative function of p, then it would be p(r, r) — > p, which would imply e(r, r) — ► 
0, that is a contradiction. 

As in section 13.11 for physical reasonableness we suppose the initial energy 
e(r, r) (and therefore p(r, r)) to be a non increasing function of r. This implies 
that we can consider, without loss of generality, the case in which also B (l) ^ 0. 
In fact, if -Bo(l) vanishes by (12.51) it has to be p(r, r) — > as r — > + . This fact, 
recalling the assumptions made on the pressure, shows that p(r, r) (that is a non 
increasing function of r) must be identically zero. But p (and therefore p) must di- 
verge at the spacetime singularity, and therefore there exists an hypersurface, such 
that p is non zero but the energy e is still regular, where we can re-assign the initial 
data on. On this hypersurface, the pressure must converge to a finite non-zero value 
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as r — > + . Then we will suppose B (l) ^ 0. Finally, we note that positivity of 
pressure on the data further implies that A (l) > and B (l) < 0. 

We are now ready to investigate lower order terms in the mass function. Recall 
that regularity of pressure along the initial data implies that ^ cannot contain first 
order terms (see (14.51) ). Then, as in (13.121) we set 

V(r,R) = V I ; ,,'•'/>'' • •••• 

i+j=3 

We now denote by e > the limit lim T _ >0 + e(r, r). By O of Remark |4~2"1 there 
exists a unique p > such that e(p ) = £q, and clearly p = hm r ^ p(r, r). We 
also denote by (3 the positive number ^(po)- Using (13.101) at r = 1 we have 

(4.8) £(r) = _|_J_ + ..., 
plus higher order terms. Observe now that, for a fixed r, 

x Bq(t) 1 /*21 *12 \ 

= + • • • = + 2 ~T + 311/03 j " Po(r) - 

If (^21, ^12) 7^ (0, 0), then po(r) — > 00 as r — > 0, and so p(r, rr) (and therefore 
j^(p(r,rr))) is sufficiently large, for r near to 0. This leads to a contradiction, 
since using (14.81) in (14.71) shows that K(r, rr) = j^u^n, —\ , for some non-zero 
constant \i independent of r, but (13.101) must hold. Then, again, B (t) = B (1)t 2 , 

andp (>) =Po(l) = Po > 0. 

Then the above argument shows that the lower order terms of the mass have the 
structure, analogue to (13.111) . 

(4.9) *(r, i?) = - ( r 3 - R 3 ) +.... 

2 V e o + Po J 

We now proceed analyzing the nature of the singularity forming at the center. With 
arguments similar to Lemma 13 .51 opportunely modified, it can be checked that 
some of the curves R x = xr 3 , for x > h sufficiently near, are emanating from 
the central singularity (if seen in comoving coordinates). Indeed, we first observe 
that, in the case of a barotropic equation of state, d2.2dl ) yields — Au = < ||, 
where the inequality is given by dominant energy condition e — p > 0. This implies 
e _I/ < y/p, and so the counterpart for the integral in (13.201) in this case has the 
following upper bound 
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Taking into account (13.221) (that still holds) in (14.51) to evaluate p(r, r 3 r), one can 
see that the integral above takes a similar form to (13.211 ): 



and so it is finite, as in Lemma [331 

Now, using (14.61) . one can compute both sides of the null radial geodesic equation 



for R = xr 3 , obtaining a similar expression to (I3.23I) . We only remark that in this 
case, since Y(r, r 3 r) cannot be infinitesimal as r goes to zero, then j"(p(r, r 3 r) is 
finite (see (14.71) ) and so is p(r, r 3 r). We can therefore conclude this section with 
the analogue of Theorem 13 .41 that is: 

Theorem 4.4. Under the hypotheses made on the equation of state in the assump- 
tion^^ for any choice of initial data for the Einstein field equations such that 

(1) the central singularity forms in a finite amount of comoving time, and 

(2) the Taylor expansion of the mass function is given by (14.9b . 

there exists solutions of (14.101) that extend back to the central singularity, which is 
therefore locally naked. 



Up today all analytical studies on naked singularities formation in collapsing 
matter of astrophysical interest (i.e. fluids) have assumed simplyfing hipotheses 
such as dust or self-similarity. 

We have shown here for the first time that among non self-similar barotropic per- 
fect fluid solutions, all those describing complete collapse for which the mass func- 
tion is regular in a neighborhood of the regular center up to singularity formation 
form naked singularities. Besides of spherical symmetry, this result is independent 
on any simplifying assumption. 

The problem of the classification of the data which leads to such singularities 
remains for future work. In particular, it is unclear if the set generating naked 
singularities is really of non-zero measure in the space of the data or not. 
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